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Abstract 

We develop a form factor approach to the study of dynamical correlation functions 
of quantum integrable models in the critical regime. As an example, we consider the 
quantum non- linear Schrodinger model. We derive long-distance/long-time asymptotic 
behavior of various two-point functions of this model. We also compute edge exponents 
and amplitudes characterizing the power-law behavior of dynamical response functions on 
the particle/hole excitation thresholds. These last results confirm predictions based on 
the non-linear Luttinger liquid method. Our results rely on a first principles derivation, 
based on the microscopic analysis of the model, without invoking, at any stage, some 
correspondence with a continuous field theory. Furthermore, our approach only makes use 
of certain general properties of the model, so that it should be applicable, with possibly 
minor modifications, to a wide class of (not necessarily integrable) gapless one dimensional 
Hamiltonians. 



1 Introduction 

Recently we developed a new method to obtain large distance asymptotic behavior of 
correlation functions in critical models from their form factor expansion in finite volume. It 
originates in a careful analysis of the form factor series for large but finite systems that be- 
come critical in the infinite size limit. In such a limit, form factors of local operators can be 
shown to scale to zero as a non-trivial power law of the system size related in particular to 
their conformal dimension ^ n the asymptotic regime (at large distances) and for large 
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(but finite) system size, the form factor series localizes to its critical part. Namely, it can be 
shown that only a subset^] of the states contributes to the leading asymptotic behavior of the 
correlation function. We called these states and their associated form factors critical. It turns 
out that the corresponding series over these critical states and form factors can be summed up 
thanks to a beautiful (purely combinatorial) multiple sum formula (see It is worth men- 
tioning that this summation formula has a deep group theoretical origin in the representation 
theory of the infinite permutation group [4l4l0l]. In the context of critical models, it gives the 
explicit correspondence between the microscopic dynamics and its critical (conformal) limit. 
Moreover, for integrable models, this result agrees with the direct asymptotic analysis of the 
correlation function performed from their multiple integral representations obtained in the 
algebraic Bethe ansatz framework ll-lf| and using rather sophisticated Riemann-Hilbert 



techniques [17|, ]1£ 

The result of this form factor series analysis is to provide, in a comparatively rather 
elementary way, the full leading asymptotic expansion of correlation functions in terms of the 
distance between local operators. Namely, for each oscillating harmonics quantized in terms of 
the Fermi momentum, it leads to the corresponding exact critical exponent together with the 
associated (non universal) amplitude. In particular, this approach succeeds to confirm, from 
a computation based on first principles, the full Luttinger liquid or conformal field theory 
(CFT) predictions 0-0 for the asymptotic behavior of correlation function in terms of 
the distance. It also provides in addition the corresponding amplitudes written in terms of 
specific and properly normalized form factors of local operators; this last information, being 
model dependent, is naturally out of reach of Luttinger liquid theory or CFT. It should 
be emphasized here that the results obtained in [l[ using this method do not rely on any 
conjectural correspondence of the model with a CFT or Luttinger liquid theory (or even its 



non-linear version see e.g. [3CH35I] and references therein). 

The aim of the present article is to apply this new method to the study of dynamical 
correlation functions and associated structure factors of critical models. As it is well known, 
the dynamical structure factors (Fourier transforms of the dynamical correlation functions) 



give the response functions under external perturbations of the system [36r|39|. As such 
they can be measured experimentally, see for examples 0, Effl]. It is therefore of 
immediate importance to design efficient methods to determine the values of these functions. 
Although the exact and analytic computation of the structure factors in integrable models 
throughout the full Brillouin zone is for the moment out of reach of our method (see however 



the method initiated in [491-4541] for a combined analytic and controlled numerics that fits 
extremely well the actual experimental datas), we will provide exact information for these 
functions in two situations of interest: the large distance and long time asymptotic behavior 
of correlation functions and the behavior of the corresponding structure factors near the edges 
of the spectrum. 

Let us be slightly more precise about the content of the present paper. To implement 
our method we consider here the ID Bose gas (or quantum non-linear Schrodinger model, 
QNLS) with delta function interaction. The QNLS model describes the evolution of quantum 
Bose fields ^(x,t), W(x,t), subject to canonical equal-time commutation relations, with 
Hamiltonian 

L 

H = j (d x & d x ^ + c^^^^ - h&V) dx. (1.1) 



lr rhis subset contains nevertheless an infinite number of states in the thermodynamic limit. 
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Here L denotes the size of the model, and we impose periodic boundary conditions. We 
restrict our study to the case of the repulsive regime c > in the presence of a positive 
chemical potential h > 0. In such a case, the Bethe ansatz provides the full description of 
the spectrum in terms of particles and holes excitations [55I 56| . Since the Hamiltonian H 
commutes with the number of particle operator, it can be diagonalized separately in every 
sector with a fixed number N of particles. In each of these sectors, the model is equivalent to 



the iV-body gas of bosons with 5-like interactions that was studied in [55l. l56l|. 

For this model we study various critical properties associated with the simplest zero- 
temperature time-dependent correlation functions in the thermodynamic limit L, N — > oo at 
N/L fixed. The latter correspond to the expectation values, in the ground state \ip g ), of 
products of the form 

W*,t)OW,Q) m <»l»ft»y.OI»>, (1 . 2) 

\ rg I Vg ) 

where 0(x, t) = e ltH 0{x, 0)e~ ttH stands either for one of the Bose fields ^>(x, t) or \E^(x, t), or 
for the local density j(x,t) = ^f<(x, t)^f(x, t). Due to the translation invariance of the model 
we set below x' = and t' = 0. 

The large-distance long-time asymptotic behavior of these two- poi nt functions (i.e. their 



behavior at t — > oo, x — > oo with x/t fixed) was studied recently in 57|, [58j via the Riemann- 
Hilbert technique (see also (59l-[69l| for related works). It was shown that, for certain correlation 
functions (in particular for the correlation functions of the densities), the leading asymptotic 
behavior cannot be inferred using a CFT-based method [23l427l| . The reason is that a cor- 
respondence with a CFT can only grasp effects taking their origin in the linear part of the 
spectrum. However, in the presence of time and beyond the x 3> t limit, there exists a 
saddle-point of the oscillating phase (see below) which is located in the non-linear part of 
the spectrum. The form factor approach developed in the present paper allows us to show 
that this saddle point gives non-trivial contributions to the asymptotic expansion of two-point 
functions. In particular, the form factor series localizes and again it can be computed using 
the same multiple sum formula as in 

The problem of the characterization of the behavior of dynamical resp onse functions on the 
particle/hole excitation thresholds was recently considered in 

[iJSS] within the non-linear 



Luttinger liquid approach in relation to the Fermi edge singularity effect, see e.g. 

0-0] • The 

form factor method used in the present paper allows us to derive the edge exponents and the 
associated amplitudes for the ID Bose gas directly from the algebraic Bethe ansatz study of 
the microscopic model, hence from a first principles based derivation. Again, in this case, one 
can show that, for the leading behavior of the structure factor at threshold, the form factor 
series localizes and that the corresponding quasi- critical sum (see Section[2] and Section [3|) 
can be computed thanks in particular to the multiple sum formula used in It appears that 
the correspondence of our results derived directly from the full form factor series with the one 
obtained through the non-linear Luttinger liquid model with one impurity relies on the rather 
non-trivial sub-additivity property of the spectrum on the particle and hole threshold lines. 
Namely, if this sub-additivity property had failed, the form factor series would have localized 
in a different way and the simple one-impurity picture used in the non-linear Luttinger liquid 
approach would be invalid. 

The article is organized as follows. We give in Section [2] a detailed description of the 
excited states of the model, which allows us to formulate more precisely the main steps of 
our method of summation over form factors. We then successively apply this approach in the 
two cases mentioned above: to the calculation of the long-time and long-distance asymptotic 
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behavior of two-point functions in Section [3l and to the computation of the edge exponents 
of dynamical response functions in Section [H Several auxiliary formulae are gathered in two 
appendices. 



2 Form factor expansion of correlation functions 



The form factor series representation of the two-point function (jl.2p is obtained by inserting, 
in between the local operators 0\x,t) and O(0, 0), a complete set of eigenstates | tj)') of the 
Hamiltonian: 



(O* (x,t) 0(0,0) 



£ 

\i>>) 

£ 

\i>') 



(^|Ot( a; ,t)|^)(^|O(0,0)|V 1) 



-itSc^+ixVc 



(y/|o(o,o)|y> g 



(2.1) 



In (|2.ip . £ ex (resp. Pex) represents the relative excitation energy (resp. momentum) of the 
excited state | tp' ) with respect to the ground state \ip g )- 

There are several advantages of considering a series such as (|2.ip . The first one is that the 
whole time and space dependence is gathered within a single phase factor. The second one is 
that there exist explicit determinant representations for the matrix elements ( t/j 1 | O(0, 0) | tp g } 
(see 57, 6^]) in finite volume. The main problem is to sum up the series (|2.ip . In this section, 
we explain our general strategy to perform such a summation in the large size limit using the 
precise description of excited states given by Bethe ansatz. 



2.1 Description of the spectrum and classes of excited states 

The form factor approach requires a precise description of the spectrum of the Hamiltonian, 
namely of the ground state and the excited states over which we sum up in (|2.ip . Within the 
algebraic Bethe ansatz framework, the ground state | tp g ) = \ ip({\}) ) of (jl.ip is parameterized 
by a set of real Bethe parameters Ai, . . . , Xn solution of the following system of logarithmic 
Bethe equations: 

N 

Lpo(\ J ) + ^e(\ J -\ k ) = 2ir(j-?^y j = l,...,N. (2.2) 

k=l 

Here the functions po(X) and 0(X) correspond respectively to the bare momentum and bare 
phase of the quasi-particles which, for the QNLS model, are given as 

p (A) = A, 0(A) = i log (^±4). (2.3) 

\ic — A/ 

The number N of Bethe parameters parameterizing the physical ground state of the model 
is fixed by the chemical potential. In the thermodynamic limit L — > oo, the ratio N/L has a 
finite limit D, and the ground state Bethe roots Xj densely fill a symmetric interval [—(7,(7] of 
the real axis (the Fermi zone) with a density function p(A). 

Excited states | ^({m^}) ) of (jl.ip are parameterized by sets of solutions /J^, . . . ,fie N , of 
logarithmic Bethe equations corresponding to different choices of integers i\ < £2 < ■ ■ ■ < £n ,: 

Lpofaj) + E 9 ®li ~ M = 27T ( £ J ~ ^Y^) > j = 1, . . . ,N'. (2.4) 

k=l 
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For the QNLS model, all roots of (|2.4[) are real, and are completely determined by the chosen 
set of integers £j (so that we have labeled them accordingly). The state parameterized by the 
solutions of (|2.4p with lj = j, j = 1,...,N', is called the ground state in the 2V'-sector (it 
differs from the physical ground state of the model if N' 7^ N) . 

Hence, the sum over all excited states ) in (|2.ip corresponds to a sum over all solutions 
{fi£ .} of the logarithmic Bethe equations (|2.4p . i.e. to a sum over all possible choices of 
integers l x < ■ ■ ■ < t N > (with N' = 2V if 0(x, t) = j(x, t), N' = N - 1 if 0(x, t) = <&(x, t) and 
2V' = 2V + 1 if 0(x, t) = \Ed(a;, t)). For a given excited state | ip') = | ^({W-}) )> the relative 
excitation momentum and energy with respect to the ground state are equal to 

N' N N' N 

^ ex = ^2po(fie a ) - J^Po(Aa), £ cx = ^2zo(fie a ) - ^eo(Aa), (2.5) 

a=l a=l a=l a=l 

where po is given by ([23]) and e (fi) = fi 2 - h. 

A convenient way of describing these excited states, in particular in the large size limit, is 
to use the language of particles and holes: an excited state with n particles/holes above the 
ground state in the 2V'-sector corresponds to a solution of (|2.4|) with lj = j except for n integers 
hi, . . . , h n G {1, . . . , N'} for which lh a = p a G Z \ {1, . . . , 2V'}. The integers h a correspond 
to 'holes' with respect to the distribution of integers for the ground state in the 2V'-sector, 
whereas the integers p a correspond to 'particles'. To each choice of such particle/hole integers 
(or quantum numbers) p a and h a one can associate, through the so-called counting function 
of the corresponding state (see Appendix [A] for details), sets of particle fl Pa and hole /2/j a 
rapidities (a = 1, . . . ,n). In the thermodynamic limit, these particle/hole rapidities tend to 
some finite values that we denote by n Pa and Hh a ■ 

There are therefore two possible ways of characterizing a given particle/hole excited state: 
either in terms of the integer quantum numbers p a and h a or in terms of the corresponding 
particle/hole rapidities. The former will be referred to as a microscopic description of the 
excited state, and the latter as a macroscopic description. The meaning of this terminology 
is the following: for size L large but finite, a finite deviation of the integer p a (resp. h a ), say, 
Pa — Pa + k with \k\ € iV, leads to a deviation of order L _1 for the corresponding rapidity 
j2p a (resp. flh a ) and deviations of order L~ 2 for all other rapidities. Therefore, although these 
two ways of characterizing the particle/hole excited states are equivalent as far as the system 
size L remains finite, one can find, in the thermodynamic limit, infinitely many microscopic 
configurations of the quantum numbers describing an excited state with fixed values ji Pa and 
fih a °f the macroscopic excitation rapidities. In particular, the thermodynamic limit of the 
rapidities of particles (resp. holes) may coincide, while the quantum numbers p a (resp. h a ) 
are always pair- wise distinct. 

We now define classes of excited states that play a role in the process of summation of the 
series (|2.1|) . We say that two excited states belong to the same class P if: 

(i) they have the same excitation momentum V e ^ and energy £ ex in the thermodynamic 
limit; 

(ii) they have the same number of particles and holes separated from the Fermi boundarie^, 
with the same rapidities [A Pa and fih a ^ n ^ e thermodynamic limit. 

2 A particle (resp. a hole) is said to be separated from the Fermi boundaries if the corresponding rapidity 
fi Pa (resp. flh a ) differs from ±q in the thermodynamic limit. Note that we say nothing here about the number 
of particles and holes which are not separated from the Fermi boundaries (except that point (i) should be 
satisfied). In particular, the total number of excitations is not fixed. 
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Since there is a one to one correspondence between excited states and form factors, the 
latter can also be gathered into classes in the sense of the above definition. 

It follows from this definition that, within a given class, the excited states (for N, L large) 
differ from each others not only by 1/L excitations around the (macroscopic) particle \x Va 
or hole /% a rapidities away from the Fermi boundaries, but also by the possible creation of 
different numbers of particles and holes on the right or left Fermi boundaries ±q. Being on 
the macroscopic level, one can in fact distinguish states of the same class only on the basis of 
their total number of particle-hole excitations on each of the Fermi boundaries. Two states 
of a given class with the same number of such particle-hole excitations can nevertheless be 
distinguished at the microscopic level by their quantum numbers p a and h a . 

In [1], we dealt with particular classes of excited states (resp. form factors) which were 
called critical states (resp. critical form factors). In such classes, all particles and holes are 
accumulated on the Fermi boundaries. In the present paper we also consider classes of excited 
states with several particles (holes) separated from the Fermi boundaries. We call such excited 
states (resp. form factors) quasi- critical states (resp. quasi- critical form factors). 



2.2 Summation process 



The main idea of our method is to use our partitioning of excited states into classes so as to 
split the summation of the form factor series (|2.ip into two steps: 

- (V|o(o,o)IV : ' 



\4>') 



—it£ cx +ixVc: 



{<!>• | O(o, o)| 



E E • 

P |V')6P 



-itScx+ixVc: 



(2.6) 



1. 



In (|2.6p we first compute the sum of form factors of a given class P, using the microscopic 
description of the excited states. We should in particular take into account, in the large- 
size representation of the form factors (see 0,[57,[H3])) the part that non-trivially depends 
on the quantum numbers. We have argued in [1] that this part is a quite universal object 
that, to a large extend, is model-independent. Due to its specific form (see Section [3] 
for details), the sum of quasi-critical form factors of a given class can be calculated in 
the large-size limit using the multiple sum formula of [l]. 

2. The second step consists in summing up the contributions from different classes. Ac- 
cording to the problem we consider (see Sections [3] and H|) , it is possible to argue that 
only certain very specific classes P will contribute, at leading order, to the critical be- 
havior we are interested in. At this stage of the calculation, it becomes enough, for the 
cases we are dealing with, to use the macroscopic description of the excited states, and 
to replace the sum over classes by integrals over the macroscopic particle/hole rapidities 
fi Pt h- Moreover, the obtained integrals are localized around very specific particle/hole 
rapidities ii v .h- 

Let us be slightly more specific. It was shown in 0, 0, [sH that particle-hole form fac- 
tors scale to zero in the thermodynamic limit. More precisely, if the system size L goes to 
infinity, they decrease with some negative power of L, which depends on the corresponding 
particle/hole rapidities: 



0(0,0) \ip g 



L- e (M) A ({^ h }\{p,h}). 



(2.7) 



IhM WW 

Here we have pointed out that the finite amplitude A depends both on the quantum numbers 
{p, h} and on the particle/hole rapidities {^ P) h}- This notation stresses some important 
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properties of the amplitude. Namely, if a particle or a hole is separated from the Fermi 
boundaries by a finite distance, then a finite change of the corresponding quantum number 
results into a variation of order L^ 1 of the amplitude A; in this case, one can say that A 
depends on the macroscopic rapidities fi Pt h rather than on the quantum numbers or that, in 
other words, it is enough to use a macroscopic description. On the contrary, for a particle or 
a hole on the Fermi boundary, (i.e. such that Jl Pa or fl] la tends to ±q in the thermodynamic 
limit), then a finite deviation p a — )• p a + k (resp. h a — > h a + k) of the corresponding microscopic 
quantum number leads to a finite variation of the amplitude; in this case the microscopic 
description should be used. This property of the amplitude becomes extremely important 
when we sum up form factors of a given class because this amounts to summing up excitations 
at the Fermi boundaries. 

Calculating the sum of quasi-critical form factors we should take into account L _1 -cor- 
rections to the excitation momentum and energy (l]. For L large but finite the latter can be 
presented as 

<p _ -p(o) , -p(i) c — £(o) i c(i) Co a) 

I ex — I ex ' 1 ex 1 °ex — <^ex ' °ex i \^-°) 

where the thermodynamic limits V^} = Vex ({^ P ,h}) and fJx = £c\* {{l-ip,h}), 

n n 

=£[KM P J-K^J], £g) = 5>0u Po ) -e(n h j), (2.9) 

a=l a=l 

are given in terms of the dressed momentum (|A.2|) and the dressed energy (|A.3j) of the particles 
and the holes at the thermodynamic limit, whereas the corrections Vex = Vex \{^p,h}\{P-> h}) 
and £ex = £ex {{l^ P ,h}\{p, h}) are of order L" 1 (see (|3.10|) . (|3.1ip for their explicit represen- 
tations). Note that V^x and Sex can be described on a macroscopic level and that they take 
the same value for all quasi-critical states in a given class. On the contrary, Vex and £e^ 
explicitly depend on quantum numbers, and therefore on the specific representative of the 
class that we consider. 

In principle, so as to be consistent with the above degree of approximation, one should 
also take care about the L~ 1 -corrections to the large-size representation (|2.7p . However, in 
the problems we deal with (long-time long-distance asymptotic behavior and calculation of 
the edge exponents), such corrections only produce subleading contributions to the results, 
and we can therefore neglect them. In particular, one can consider the power L~ 9 ^p< h ^ as a 
common prefactor L~® p for a given class P of quasi-critical form factors. 

Therefore, restricting the form factor series (j2.6|) to the sum over quasi-critical states 
belonging to a given class P, we obtain, at our order of approximation in L^ 1 , 



E 

IV>')€P 



(V' I o(o,o) iv 2 



L -e Pe -it£^+iM<? e^-(fe' l >)- itf -(fe' l >)^({ y u pA }|{p^}). (2.10) 

MdM 



Here the values of the rapidities //„ ^ in the arguments of the amplitude A are fixed. 

In [l|], we dealt with the sum ()2.10p in the case £^ = and t = 0. We have in particular 
shown there that, in this case, the computation of the sum eventually produces a factor L® F 
precisely compensating the prefactor L _6lp , thus ensuring that the final result has a finite 
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thermodynamic limit. It means that the sum over excitations on the Fermi surface provides 
an effective dressing of the original form factors. 

In our present case, with more general classes such that Sex / 0, this dressing compensates 
the vanishing pre-factor L~ 6p only partly. However, the remaining negative power of L 
disappears when we sum up the contributions of different classes: it can actually be absorbed 
into the integration measure when one replaces the summation over classes in (|2.6p by an 
integration over the particle/hole rapidities separated from the Fermi boundaries. 

To achieve this second step, and to evaluate the two-point correlation functions, it finally 
remains to compute these integrals over {f-iph}- Whereas their explicit analytic evaluation 
is hardly possible in the general case (this is due to the very complicated dependence of the 
amplitude A on the rapidities {[i p ,h}, see e.g. [571]), the situation can be simplified if, for 
some reason, these integrals are localized around one point (or several points): the amplitude 
can then be replaced by a constant, namely by its value at the point of localization. In the 
examples discussed below, we are precisely dealing with such a case. 



3 Long-time and large-distance asymptotic behavior of the 
two-point functions 

In this section we study the long-time and large-distance (x — > oo, t — > oo, with x/t being 
kept constant) asymptotic behavior of the time-dependent two-point correlation functions, 
using the form factor summation process sketched in the previous section. 



3.1 General scheme 

Following the same method as for the time-independent case we first need to identify, 
within the series (|2,ip . the relevant form factors contributing to the leading asymptotic of a 
given oscillating harmonics. 

In the thermodynamic limit, for states with n particles and n holes, the oscillating phase 
can be presented in the form: 

n 

-itS® + ixPj® = -it Wp«) - > 

a=l 

where 

u(A) = e(\) - !„(A). (3.2) 

In the limit x — > oo, t — > oo with x/t = const, the phase factor becomes rapidly oscillating. 
Similarly to what happens for oscillating integrals, we therefore expect that the main con- 
tributions to the sum come either from the boundaries of the summation interval (the Fermi 
boundaries ±q) or from the phase saddle-point. In the equal-time case studied in [if, there was 
no phase saddle-point, and the leading asymptotic contribution was solely issued from states 
characterized by particle and hole excitations close to the Fermi boundaries. However, in the 
present dynamical case, the saddle-point contributions should also be taken into account. 

The existence and the uniqueness of the saddle-point of the phase u(X) (|3.2p is closely 
related to the properties of the sound velocity v(X) in the Bose gas: 

de e'(X) 

- w P - M (3 - 3) 
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It is believed that this function is strictly monotonic: v'(X) > 0. This property can be proved 
if the Hamiltonian parameters c and h are such that c/q(c, h) <l 0.6. Otherwise, to the best 
of our knowledge, the corresponding proof is missing up to now (see e.g. [Zf3]). The numerical 
analysis of equations ()A.2p and ()A.3p does nevertheless confirm that v'(X) > for arbitrary 
values of c and h, and we therefore assume this property. It also follows from equations (|A.2j) 
and (|A.3p that v(X) — > ±oo as A — > ±oo. This ensures that the phase u(X) has a unique 
saddle point Ao determined by the equation v(Xo) = x/t. Moreover, the monotonicity of v(X) 
associated to the fact that p'(X) > (c/ (|A.2|) ) also ensures that u"(Xq) > 0. 

Introducing the Fermi velocity v F = v(q) = —v(—q), we can define two asymptotic 
regimes^: 

• the space-like regime (\x/t\ > v F ) which corresponds to a saddle-point outside of the 
Fermi zone |Ao| > q; 

• the time-like regime (\x/t\ < v F ) which corresponds to a saddle-point inside of the Fermi 
zone |Ao| < q. 

We now describe the excited states (and the corresponding form factors) giving rise to 
the leading asymptotic terms of the correlation functions. In the space-like (resp. time-like) 
regime, the main contribution obviously comes from form factors associated to states for which 
the particles (resp. holes) which are separated form the Fermi boundaries (if any) are close 
to the saddle-point Ao- In other words, a contributing state with n particle-hole excitations 
in the space-like (resp. time- like) regime is such that there are no particles (resp. holes) with 
rapidities {^ Sa } lying in a small vicinity of the saddle point Ao, the other n — no particles (resp. 
holes), as well as the total number n of holes (resp. particles) being on the Fermi boundaries 
±q, with no > 0. To describe such states, it is convenient to introduce a parameter r such 
that r = 1 in the space-like regime and r = — 1 in the time-like regime. Let also n^ be 
the number of particles with rapidity ±q, and be the number of holes with rapidities ±q. 
Then, in both regimes, the above-introduced numbers counting the amount of particles (resp. 
holes) all sum up to the total number of excitations: 

n+ + n~ + - ^ T ^ n = n+ + n~ + - - T ^ n = n. (3.4) 

Let I refer to the right particle/hole discrepancy number t = n+ — ri£. Due to the constraints 
(|3.4p . this means that the left particle/hole discrepancy number can be written as Tip n^ — 
— TUq — i. 

The states described above form a class of what we call quasi- critical states. This class is 
characterized by the integer numbers £, r and no, and by the rapidities fi si , . . . , fj, SnQ localized 
in a small vicinity of Ao with no > 0. If these parameters are fixed, then, in the thermodynamic 
limit, all states within this class have the same energy and momentum, namely 

n n 

^=rJ2e(p 8a ), V^=rY / p^ Sa ) + (2£ + rn )k p , (3.5) 

0=1 a=l 

where k F = p(q) is the Fermi momentum. We refer to this class of quasi-critical states (and 
form factors) as the P T n / class. 

As mentioned in Section [21 all the form factors of the P rno i class scale to zero with the 
same exponent that we denote by Tno ,g- The latter was computed for generic particle-hole 

3 The analysis of the asymptotic behavior for x/t = ±v F goes way beyond the scope of this paper. 
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form factors in 0, It was shown there that it can be expressed in terms of the values of 
the shift function (see (|A.13[) ) on the boundaries of the Fermi zone: 



with 



-rn , 



F Tno/ (q) + AN. 



(3.6) 



(3.7) 



The integer AiV = N' — N is for the correlation function of densities and ±1 for the 
correlation functions of fields, and the shift function F Tno £ associated with the P rn0i £ class 
takes the form 



F Tno ,eW 



-AN 



Z{\) 



no 



+ cf>(\,q) -^(A,g) + (£ + mo)^(A,-g)-rJ^^(A,MO (3i 



o=l 



in terms of the dressed charge Z{X) and dressed phase <j)(\,v) defined respectively in (|A.5[) 
and (jXi]). 

In order to fix a quasi-critical state inside the P rno t class, we need to introduce the 
quantum numbers labeling the particles and holes with rapidities equal to ±q. We use the 
standard (see [l[) re-parametrization of the original quantum numbers {p} and {h}: 



Pj =Pj+N', if n P 

Pj = l-Pj, if Hp 

hj = N' + l-hf, if Hh 

hj = hT, if fi h 



(3.9) 



These integers allow us to express the amplitude A and the finite-size corrections to the 
oscillating phase. 

A careful analysis of the finite-size corrections in (|2,5p enables us to show that the correc- 
tions to the momentum and energy in (j2.8|) have the form 



V^i — —V -I 

'ex ~~ £ ' TTl ,£ I 



E (p* - x ) + E - XX - 1) - E h « 



0=1 



a=l 



rno, 



+ 



a=l 
n7> 



o=l 



E^ - 1) + E ^ + E^~ - !) + E ^ 



0=1 



0=1 



0=1 



a=l 



(3.10) 



(3.11) 



where V Tno> £ and £ Tn0: e are constants (at least at the leading order in L) for all states belonging 
to a given P T no,£ class. 

The explicit form of the amplitude A for the quasi-critical form factors can be obtained 
following the same lines as for the critical form factors. There is only one subtlety due to the 
localization of the rapidities ^ Sa in the vicinity of Ao- The matter is that, due to the fermionic 
structure of the excitations, the form factors vanish as soon as two or more rapidities of 
particles coincide. Actually, the amplitude A has the form A = A 2 ({fi Sa })A, where A({fi Sa }) 
is the Vandermonde determinant of the rapidities {n Sa }- The effective amplitude A has a 
non- vanishing limit as fx Sa — > Ao- Thus, as a first approximation, we can set all {n Sa } equal 
to Ao in the arguments of A, but we should keep the prefactor A 2 ({fi Sa }) as it is. 

So as to express this amplitude, it is convenient, as in the time-independent case, to 
introduce the simplest form factor of the ~P rn0t e class. The corresponding quasi-critical state 
I V'rno t ) * s defined as follows: 
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• the no rapidities fi Sa (which are the same for all representatives of the class) are in a 
small vicinity J\ of the saddle point Ao; 

• the distribution of particles and holes on the right Fermi boundary is as follows: if £ > 0, 
there is no hole (ri^ = 0) and there are £ particles (rip = £) which are characterized by 
the integers p+ = a, a = 1, . . . , £; if £ < 0, there is no particle (n+ = 0) and there are \£\ 
holes (nt = —£) which are characterized by the integers h£ = a, a = 1, . . . , —£; 

• the distribution of particles and holes on the left Fermi boundary is as follows: if £+tuq < 
0, there is no hole (n^ = 0) and there are —I — tuq particles (n~ = — £ — too) with 
p~ = a, a = 1, . . . , — £ — tuq; if £ + tuq > 0, there is no particle (n~ = 0) and there are 
£ + too holes (n^ = £ + too) labelled by the integers h~ = a, a = 1, . . . ,£ + too- 

Using this representative of the P T „ 0j £ class, we introduce the renormalized form factor 

LF^ = hm hm L" Tn a- 1 A (\a s }) , "\ , — ; ; — ; -. (3.12) 

Then, for large L and after setting fi Sa = Ao, the leading order of A can be written in the 
following form: 

I{moA/r\ \\[„± h ±\\-\ T ,2 ° 2 ( l + F rn ,e) G ' 2 ( 1 ~ F rn /) 

* R n},„ t ({P + }dh + }\FX,J R„-, n -({p-}.{h-}\ - F T - n J , (3.13) 
where G(z) is the Barnes function satisfying G(z + 1) = T(z)G(z), and the function R n ^ n i is 



n n 

/ ■ T-i \2n' 

r a \ sh\\r?\ ( sm7TF \ i>t TT r Pfc + F ) tt r fofc ~ F > 



i=i fc=i 

Thus, the sum of the form factors of the P Tno l class reads 



-it£ e x+ixVc: 

l*'>6 p rn ,f 



(V' 1 0(0,0)1^ 



2 



x ^ e 4Xp - ) - it£ -i( rno ^({Ao}|{p ± ,/i ± }). (3.14) 

Using the explicit representations ()3.10p . ()3.11|) for and and (|3.13p for A^ Tn °'^ we 
see that the sum over the quantum numbers {p 1 * 1 }, {h^} can be factorized into a sum over the 
integers {p + },{h + } on the one hand, and a sum over the integers {p~},{h~} on the other 
hand. Both of these sums can be computed explicitly via the identity used in [l| (see also 

/ 

n n 

EE E W ^ ^ Rn,n'({Ph{h}\ F ) 

n,n'>0 l<pi< --<p n <oo l<hi< - <h /<oo 



n—n =r 



G 2 (l + r + F) yf^-m 

G 2 {l + F) (l- W )(F+ry [SAb) 
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In our case one should set 



rno, 



for the sum over {h + }, and 



™ = e -¥(*+M) 



rno, 



rn , 



(3.16) 



(3.17) 



for the sum over {p~}, {h~}. Moreover, since the identity (|3.15p is valid for < 1, we should 
regularize the form factor series by adding a negative imaginary part to the variable t, i.e. by 
replacing t by t — iO. Note that such a regularization is very natural since it increases the 
convergence of the original form factor series f|2. 1 1) . Then, the sum over form factors of the 
P T „ ^ class results into 



E 



-itE ex +ixVe: 



IV>'}eP rno ,f 



< V' I 0(0,0) ] v s 



\T° \ 2 A 2 ({u s \)L- n o e ixV ™- it£ ™ 



X e X^-^F'K^- 1 )-x( X +^F*)^+™o)(^+™o + l)+i;(^ no/ -^ Tno ,f) 

It is quite remarkable that now the thermodynamic limit can be taken in the last two lines 
of this equation. Obviously the limit of the second line is 1. Using the explicit form of the 
exponent Tn(h e (|3.6[) one can easily see that the last line has a finite thermodynamic limit: 



lim L~ e ™0' l+no 
L— >oo 



1 _ e ¥(-^*))" (F ^' f+ ' )2 (l _ e -^(x + , Fi )^^o^ + ™o) : 

g if [sgn(x- ?)F i)(F+ no ^+£) 2 -sgn( :C + ?)F i)(F T - noi ,+£+rno) 2 ] 

" |2tt(x - v F t)\ {F ™^ +l)2 \2it(x + v F t)f^ +l+Tno)2 



. (3.19) 



The first line in (|3.18|) is proportional to L~ n °, but this factor will be absorbed into the 
integration measure in the second step of the calculation, when we sum up the contributions 
from different P Tno / classes. Therefore in the large L regime, we can safely take the limit 
L —7- oo in the last two lines of (|3.18p . which gives us 



E 

l-/>')€P rn0i 



-it£ cx +ixV c: 



o(o,o) |W) 



O |2 a2 



A 2 ({fi Sa })L- no e lxV ™ - lt£ ^ 



if [sgn(x-v F t)(F+ no te +e) 2 -sgn(x+v F t)(F- noe +e+Tn ) 2 ] 

X t ; r. (3.20) 

|2vr(z - v F t)\^o^ 2 \2n(x + v F t)f™o^™^ 

Thus, we have calculated the sum over quasi-critical form factors of a single P T n / class. As 
we have seen, most of the singularity (apart from the L~ n ° factor) of the form factors has 
been absorbed by this process. In fact, such a summation can be interpreted as a dressing 
of the original bare simplest form factor of the class by a cloud of particle-hole excitations 
around the Fermi zone. The resulting dressed form factor has now a macroscopic description 
in terms of the rapidities [i Sa of the particles and holes separated from the Fermi boundaries. 

The second step of our computation consists in summing up such contributions from 
different classes. This means that we should sum up with respect to all possible values of the 
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non- negative integer no, of the integer £, and also of the rapidities fj, Sa localized in a small 
neighborhood J\ of the saddle point Ao- The latter sum can be replaced by an integral in 
the L — >■ oo limit according to the rule: 

T /(/O T — > [ /(/ x )p(/ u ) d / u for any regular function /, (3.21) 

L L — ' L->+oo J 

where p(/i) is the density (jA.ip . Since we have a multiple sum over the tiq variables fj, Sa , 
we obtain an no-fold integral and, as we claimed before, the prefactor L~ n ° in (|3.20p does 
disappear: 



lim ^- n °A 2 ({n Sa })e 



n l 



/no 
A 2 (M) e^-^ H p ^ a ) dfx a , (3.22) 
a=l 



where vix and are given by (|3.5|) . Using standard saddle-point considerations, we expand 
the oscillating phase up to the second order term around Aq: 



ixpW - it£<® = ix(2e + rn )k F + irn (xp(X Q ) - 1 e(A )) - itr^^- - A ) 2 , (3.23) 



0=1 



and, using the fact that u"(Ao) > 0, we reduce our integral, in the t — > oo, x — > oo limit, to 
the Gaudin-Mehta integral: 



1 p n o ii 

W- J ■ LJ r 

t a=l 
->a 



p no (A )e- ir f n o s g n W /■ 2 -ft |t|«"(A n ) 2 



n ! 



/A^ne-^^^l + Olt- 1 )) 

m 0=1 



|to"(Ao)| n o/ 2 V ; 



Combining (|3.20p with ()3.24p and taking the sum over no and £, we obtain the leading 
asymptotic terms for all the oscillating harmonics of the two-point function at x,t —> oo, 
x/t = v(X ): 



oo oo 



n^O^col^'W-^'Wro/ 2 

^^(no.£) '^oJ 2 exp N^+^o)^ +irn (xp(X ) -te(X ))] 



where 



n 2 



V? T (n , *) = sgn(x - + £ ) " sgn ( x + v F t )( F rn ,i + i + Tn o) ~ r sgn(t)-^. (3.26) 
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We recall that r = 1 in the space-like regime and r = —1 in the time-like regime, and that 
the values F Tno t of the shift function are given by (|3.7|) . (|3.8|) . where one should set fj, Sa = Xq. 
Using the fact that \x/t\ > v F in the space- like regime and that \x/t\ < v F in the time-like 
regime, one can slightly simplify (|3,26l) : 

<p+(n J) = sgn(x) [(F+ o/ + if - (F~ o/ + i + n f] - sgn(t)^, (3.27) 

2 

p_(no,*) = -sgn(i) [(F+ n0j£ + if + {FZ na + i - n f - ^] . (3.28) 

We keep in our result the sum over all possible values of uq and I because of two reasons. 
First, each of these terms is associated with a different oscillating harmonic (and is leading in 
the corresponding harmonic). Second, generically it is not easy to determine the leading term 
in the series (|3.25|) . It may depend not only on the specific operator O and the position of 
the saddle point Ao, but also on the original parameters of the model (the coupling constant 
c and the chemical potential h). 

In the remaining part of this section, we specialize this general result to the case of the 
field-conjugated field and of the density-density dynamical correlation functions. 

3.2 Correlation functions of fields 

The asymptotic behavior of the correlation functions ( t) ^(0, 0) ) and ( W(x, t) ^(0, 0) ) 
can be obtained directly from the general scheme presented in the last subsection. 

The form factor sum for the two-point function ( ^(x, t) W (0, 0) ) involves quasi-critical 
form factors with AN = 1 . Using ()A.13p , ()A.6|) , (|A.7|) , we can rewrite the exponents of f|3.25|) 
as 

F rn ,e + e = £Z+^+rn ^ + , (3.29) 

F ™ ,e + i + Tn o = iz -Jz+ Tn o®-i ( 3 - 3 °) 
where Z is defined as the value of the dressed (fractional) charge Z(X) on the Fermi boundary: 

Z = Z(q) = Z(-q). (3.31) 
We have also introduced the following notations: 

$+ = <f>(q, -q) - A ) = -<f>(q, Ao) + \ (3.32) 

<D_ ee 1 + 0(-g, -q) - 4>{-q, A ) = A ) + - + ^ - , (3.33) 

where <fi(\, fj) is the dressed phase (|A.4|) . and we used the identities ()A.6p and (|A.7|) 

Now the general formula can be directly applied to give the following leading (in each 
harmonics) asymptotic behavior, 



(*(» l t)*t(o,o)> = £; f: (^p(Ao)r^(i + n 



n =o i=-oo i te"(X )-x P "(X )\^ 
2n{x-v F t)f z+ ^z +™o *+) 2 |2vr(x + w F t)|^-2^ +Tn ° * - ^ 



is<Pr(no i) Sll! exp [ ix ( 2i + Tn ^ k F + irn (xp(X ) -te(X ))] 
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To get some insight into this asymptotic behavior, let us look at the possible values of the 
critical exponents S n0) i jT) \ associated to the terms behaving as |x| -<5n o^ T ' A o in (|3,34|) . They 
are positive numbers given as the sum of three squares 

<W,Ao = Y + ( £Z+ 1Z + Tn ^+f + { eZ ~JE + rn o$-) 2 . (3.35) 
For no = 0, the critical exponents do not depend on r, Ao and we get 

6 0>i = 2£ 2 Z 2 + ^. (3.36) 

Hence, the no = series over £ in our result (|3.34p fully reproduces the Luttinger liquid 
theory and CFT predictions and provides in addition the exact values of the associated am- 
plitudes. In this case, the dominant term is given obviously by £ = with critical exponent 
1/2Z 2 < 1/2 as in the QNLS model with finite positive coupling Z > 1 {Z = 1 only in the 
free fermion limit c — > +oo). 

The rest of the series, namely the terms with no > 1, corresponds to the contribution of 
the saddle point located in the non-linear part of the spectrum (separated from the Fermi 
boundary). It is important to stress that CFT and Luttinger liquid theory fail to predict 
such contributions and the corresponding oscillatory harmonics which explicitly depend on 
the saddle point Ao (hidden in t<J>±). 

To see the possible relevance of these additional terms, let us discuss now the relative 
values of the different critical exponents. For no > 1, S no ^ T \ Q > 1/2, whereas 1/2Z 2 < 1/2 
(it is equal to 1/2 only at free fermion point). Hence, the term no = 0, £ = gives the 
dominant contribution to the asymptotic series (|3.34p . which confirms the CFT prediction. 

The origin of the first subleading term of the series ()3.34p is more involved. In the CFT 
and Luttinger liquid theory it is given by no = 0, £ = ±1 and the corresponding critical 
exponent 2Z 2 + 1/2Z 2 is greater than or equal to 5/2. It so happens that the leading saddle 
point contribution can be dominant with respect to this first subleading CFT contribution 
(see also [58j). 

To see this, let us consider a very simple example, namely the free fermion point where 
Z = <£_ = 1 and <!>+ = 0. In this case, all the dependence of the critical exponents fino,£,r,\o 
on Ao reduces to the value of r = ±, namely it just indicates whether we are in the space-like 
(r = 1) or in the time-like (r = —1) asymptotic regimes. We have for the free fermion point 
critical exponents, 

SnoAr = Y + ( i+ \T + ( i ~\ + Tn °) 2 - (3 - 37) 

For a given no, $n ,i,T as a function of £ (considered as a real variable) would reach its minimal 
value at £ = — tuq/2. Hence for no = 1, the dominant term is realized with the nearest integer 
values £ = 0, — r with corresponding critical exponents <$i,q,t and 5i- T>T . 

In the space-like regime, r = 1, and the leading critical exponents coming from the saddle 
point are <5i,o,i = 1 and <5i,-i,i = 1, which produce two dominant terms in (|3,34p with respect 
to the first subleading CFT critical exponents #o,±i,i = 5/2. Hence in this case, the saddle 
point provides the first oscillatory subleading contribution to the asymptotic behavior of the 
series (|3.34p . 

In the time-like regime, r = —1, and the leading critical exponents coming from the 
saddle point are <5i,o,-i = 3 and = 3. Hence, they are subdominant with respect to 

*o,±i,-i = 5/2. 
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The two-point function ( ^(x, t) \I/(0, 0) ) can be treated in a very similar way. The rele- 
vant form factors are those with AN = — 1. The corresponding exponents in this case are 

Frn ,i + £ = tZ-7^+™ <S> + , (3.38) 
Knoi + Z + ^o = «+l + ™ $-, (3.39) 
leading to the following asymptotic expansion of the correlation function: 

(V2^p(\o)) no G(l + no) 



t)*(0,0)) = ^ 

no=0 £=— oo 



oo 

^ 't£"(A ) -xp"(A )| n o/ 2 



i\ <p T (n ,l) ^rj 2 ex P [ ix ( 2£ + Tn o) k F + ^n (xff(A ) - t e(A ))] 

' |27r(a;-i; Ji ,t)|^-^ +Tn o*+) 2 |27r(x + w Ji ,t)|^ + ^ +Tno *-) 2 ' 

the phase ip T (no,£) being now built out of the shift functions (|3.38|) . ()3.39p . 

The analysis of this asymptotic series is similar to the one given above for f|3.34|) . The 
CFT and Luttinger liquid predictions are reproduced for no = and the term with no = 
and I = is dominant with critical exponent 1/2Z 2 < 1/2. As for the first subleading term, 
the discussion is very similar to the previous case with one notable difference: taking the 
example of the free fermion point, the first subleading oscillatory contribution is provided by 
the saddle point, with no = 1, in the time-like regime, while in the space-like regime it comes 
from the terms no = and £ = dbl. 

3.3 Correlation function of densities 

The asymptotic behavior of the time-dependent density-density correlation function can also 
be inferred from the general framework developed in Section 13.11 There is however some 
peculiarity in this case due to the fact that F e = for £ = no = 0. As a consequence, the 
method described above, when applied directly to the two-point function, only enables one 
to obtain the trivial, constant part D 2 of the non-oscillating contribution to the asymptotic 
expansion (the oscillating terms do not suffer from such effects). This apparent problem can 
be circumvented by considering, just as in the time-independent case, a small modification of 
the initial form factor series. 

Namely, as in [l| (see also (stJ), it is convenient to introduce the generating function 
( e ltH e 2ma Q( x ) e ~ ttH a ^ w here a is a real number (twist). Here we define the operator Q(x) as 

X 

Q{x) = J j(y)dy, (3.41) 
o 

and H a corresponds to the quantum non-linear Schrodinger Hamiltonian (jl.ip for which we 
impose twisted quasi-periodic boundary conditions ^{x + L) = e 2ma ^{x) instead of peri- 
odic ones. The time dependent correlation function of densities can be obtained from the 
derivatives of this generating function (see Appendix |B|) : 



(j(x,t)M0)) = - ^^ 2 (e UH e 2 -^e-^ 



(3.42) 

a=0 



The form factor approach described in Subsection 13.11 can be very simply adapted to 
the evaluation of this generating function by introducing the eigenstates \ip' a ({/*})) of the 
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Hamiltonian H a with quasi-periodic boundary conditions. These eigenstates are given by the 
solutions of the twisted Bethe equations: 

N 

L Po @e j ) + ^9(fi e .--ii ik ) = 2ir(l j + a j = l,...,N'. (3.43) 

k=l 

One can show (see Appendix [B]) that the generating function has the following expansion over 
the complete set of twisted eigenstates | ip' a ({fi}) )'■ 

( e itH e 2niaQ(x) e -itH a ^ = e - ^^ {J* ({$})) \ (3-44) 

^(M)> 

where 

{m) H' a m)\\u 9 \\ mmnwM' { ' ; 

and £ ex and Vex are defined as in (|2.5|) in terms of the rapidities ju^ o solution to (|3.43|) . 

The general scheme of Section [2] applies to the computation of the asymptotic expansion 
in the a-twisted case with only slight modifications. The shift functions F^ nQ £ can be directly 
obtained from the Bethe equations. The corresponding exponents are 

F+ m/ + £=(a + £)Z + rn ^ + , (3.46) 
F rn ,i + £ + Tn = {a + l)Z + rn (3.47) 

and the phases (p T (riQ,£) take the following form according to whether we are in the space-like 
or time- like regime: 

2 

ip+(n ,£) = sgn(x)[((£ + a)Z + n r$ + ) 2 - {{£ + a)Z + n r$_) 2 ] - sgn(i)^, (3.48) 

2 

ip-(n ,£) = -sgn(t) [((£ + a)Z + n r$+) 2 + {{£ + a)Z + n r$_) 2 - ^] . (3.49) 
It leads to the following asymptotic behavior for the generating function: 

, e itH e 2niaQ(x) e -itH a)= y^ (V2^p(X )) no G(l + ftp) 

n^O^I^"(^)-^"(Ao)ro/ 2 

I -pa 1 2 ix(2£+2ct+nQT)k F +inQT (xp(Xo)—t e(Xo))+i^^r(no,£) 
I rriQ ,^ I 

X (2vr|x - t; F t|) ((£+a)2+nor * +)2 (2vr|x + ^|) («+°)*+«o-*-) 2 



(3.50) 



Here J-^ is the thermodynamic limit of the properly normalized (cf (|3.12p ) a-twisted scalar 
product (|3.45p for the specific solution of ()3.43p which corresponds to a choice of the integers 
£j as described on page [TUJ 

It then remains to differentiate (|3.50p with respect to a and x so as to get the asymptotic 
behavior of the two-point function. Just as in the time-independent case, one should consider 
the case £ = no = separately. Indeed, \J- a | 2 has a second order zero at a = except in 
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the case I = uq = for which |J^ | 2 = 1 at a = (see [l] for details). Using (TRlH) we obtain 

( jix two o)) = x2 + vlt2 +T T ( ^ (Ao))n ° G(1 + no) 

\t\+n >0 

\P ^ e ix(2e+nor)k F +iTno(x p(X ( ))— t e(Ao))+if</3 r (no,f) 

X (2 7 r|I-, F t|) ( ^ + ™ 0$+)2 (2 7 r|x + , F t|)^ + ™°*- )2 ' ^ 

We stress that, on the level of (|3.5ip one should set a = in the expressions (|3.48p . ()3.49p 

for (p T (n ,£). 

Let us look now at the possible values of the critical exponents 5n ^,r,A associated to the 
terms behaving as |:c| — * n o>^' r > A a in (|3.51|) . Besides the constant term, the non-oscillating term 
has critical exponent equal to 2. For oscillating terms in (|3.5ip the critical exponents are 
given as the sum of three squares 

ng / n2 / ^2 



<W,r,Ao =2 + V Z + Tn °®+) + \ tZ + Tn ^~) ■ ( 3 - 52 ) 
For no = 0, they do not depend on r, Ao and we get 

5 0/ = 2£ 2 Z 2 . (3.53) 

Hence, as for the correlation functions of the fields, the series over I obtained for no = 
reproduces the CFT and Luttinger liquid predictions, with the exact values of the amplitudes 
given here in terms of specific (renormalized) form factors of the density. The dominant term 
of the asymptotic expansion (|3.5ip is evidently the constant term. The origin of the first 
subleading term is more involved in general. 

The terms with uq > 1 are subdominant with respect to the non-oscillating contribution 
which has critical exponent equal to 2. Hence, after the constant term, the leading behavior 
is produced either by the uq = and hq = 1 terms or by the non-oscillating term. For 
no = 0, the leading terms correspond to t = ±1 with critical exponent 2Z 2 > 2, hence greater 
than or equal to the critical exponent of the non-oscillating term, the terms with no = and 
I > 1 being evidently subleading. It is more complicated to select the leading contribution 
for no = 1 and to compare it in general to 2 or 2Z 2 since it depends on the regime and on 
the values of the constants $±. 

To see at least the possible relevance of these different terms, let us consider again the 
example of the free fermion point. There, the leading non-constant terms correspond to 
no = 1, I = and uq = 1, £ = — r, both with the critical exponent 3/2, which is dominant 
with respect to the non-oscillating term and to the first subleading oscillating term predicted 
from CFT (both with critical exponent 2). In this example, the saddle point provides the 
main contribution after the constant term a soon as x/t takes a finite value in the asymptotic 
regime. Hence, in this case, Luttinger liquid theory and CFT fail to predict the leading 



asymptotic behavior of the density correlation function (see also 57]). 



4 Behavior of dynamical response functions on the particle/hole 
excitation thresholds 

In this section, we use the form factor expansion to study the dynamical zero-temperature den- 
sity structure factor S(k,co) (DSF) and spectral function A(k,u), which are defined in terms 
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of the space and time Fourier transform of the two-point dynamical correlation functions: 



S(k,u)= / e i(ult ' k ^ {j(x,t)j(0,0))dxdt, 



(4.1) 



A(k,u) 



— sign(w) R f e^"*-**) ( T*(z, i) *t( 0j o) ) dx dt. 



(4.2) 



E 2 



Here T denotes the time-ordering operator. S(k,tu) describes the probability to excite the 
ground state with momentum and energy transfer (k,co), whereas A(k,u) describes the re- 
sponse of the system to the addition of a physical particle or hole with momentum k and 
energy to. Both quantities can be measured experimentally 



4.1 General scheme 

Let us consider the space and time Fourier transform of a two-point function of the type (jl.2p . 
When evaluating this quantity at some point (k,uj) from its form factor series, we see that the 
only contributions to the form factor sum come from states | ip' ) with excitation momentum 
and energy precisely given by (k,uj): 

J dxdte^-^i O^i) 0(0,0)) = ~ £ex) S(k - V ex ) 

This is due to the fact that, in the form factor sum (|2.ip . the only x and t dependence is 
contained in the phase factor. Note in particular that the quantity f|4.3[) vanishes for uj < 0, 
becaus^fl £ cx > 0. It is also invariant under the change k — > — k, since the existence of a state 
with Vex implies the existence of a state with — "Pex- This means that we can restrict our 
study to the domain k > 0, lo > 0. 

We shall consider the behavior of the above quantity in the vicinities of Lieb's one-particle 
or one-hole excitation spectra [if]]. In other words, let A denote the rapidity of an arbitrary 
particle (A > q) or an arbitrary hole (A e] — q, q[) located at finite distance from the endpoints 
±g of the Fermi zone. Let k p = p(\) — k F (resp. k^ = k F — p(X)) be the momentum of the 
excitation corresponding to a particle at A and a hole at q (resp. to a hole at A and a particle 
at q), the corresponding energy being equal to e p = e(A) (resp. to = —e{\)). We are 
interested in the 5u — > behavior of fj4.3|) when 

k = k p = p(X) — k F and oj = e v + 5ui = e(A) + 5uj (for A > q), (4.4) 

i.e. in the vicinities of the one-particle excitation threshold, or when 

k = kh = k F — p(X) and lo = + 5u = — e(A) + 5uj (for A e] — q, q[), (4.5) 

i.e. in the vicinities of the one-hole excitation threshold. 

Recall that we assume the sound velocity v(X) to be a strictly monotonic function: v'(X) > 
0. Then one can easily see fldt\ that a given combination of momentum and energy (k p ,e p ) 
(resp. (kh, Eh)) belonging to the threshold can only be realized by a unique particle/hole con- 
figuration at the macroscopic level: one particle (resp. hole) with rapidity A, one hole (resp. 

4 Note however that, due to the presence of the time-ordering operator, the spectral function 1)4.211 does not 
vanish for uj < 0. 



(y/|Q(o,o)|V> g 
IIV-'ll IhM 



(4.3) 
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particle) on the right Fermi boundary +q. On the microscopic level, however, one can dress 
this single state by an arbitrary number of additional particle-hole excitations with rapidities 
accumulating on the two endpoints ±q of the Fermi zone, so that the total momentum and 
total energy of these additional particle-hole excitations is effectively zero at the thermody- 
namic limit. It means that we should consider quasi-critical states of the P±i=pi class, with 
only one rapidity fj, separated from the Fermi boundaries. If we admit a small deviation of 
the energy e p + Soj (resp. Eh + 5uj), then p belongs to a small neighborhood J\ of A, tending 
to A in the 5uj — > limit. Hence, we should sum up the contributions of the corresponding 
P±l,=pi classes, for all such values of p. 

In fact, all this was already done in Section [3] for a more general case. The only difference 
is that now the ratio x/t is not fixed and the point A is no longer a saddle point. Therefore, 
to take into account the variation of the oscillating phase in the small neighborhood J\ of 
A, it is no longer necessary to develop the function u(p) up to the second order as in (|3,23p . 
Instead we simply linearize it around A: 

xp(fx) - te(fi) = xp{\) - te{\) + [xp'(X) - te'(X)] - A) + . . . . (4.6) 

In all other respects, one can compute the contribution of the P±i,=pi classes form factors by 
the method described in Section [3) 

For example, the contribution of the form factor sum in the vicinities of the particle 
threshold (that is the contribution of the sum over the corresponding Pi _i classes) is 



(OHx,t)O(0,0)) kp>Ep+SuJ = lim V L-W 

L— i>oo * — ' 



(A) Gu-A)-ite'(A) (m-A) 
£-01,-1+1 \ jrO I 2 e -ixk F +ixp(X)-ite(X) 



, (4.7) 



[i _ e ¥(-M)] (F it-i- 1)2 [l _ e -¥(-+M)] (F w )2 

Recall that, as discussed previously, one should use the regularization t — >• t— iO when applying 
the summation identity (|3.15p . Taking the thermodynamic limit in (|4.7p . we obtain 

(OHx,t)O(0,0)) kp , £p+Sw = P (X) I d Me ^W(A^A)-<te'(A)^-A) 

I -pO I 2 p — ixk F +ixp(\)— ite(\) 
r 1,-1 1 



•/a 



[-2m(x - v F t)] {F ^~ 1)2 [2vri(x + v F t)] {F ^ )2 ' ^ 

where we have used (|3.2ip and replaced p(//) by its value p(X) at the localization point. In 
contrast to (|3.19p . we prefer to present the result in terms of x ±v F t instead of their absolute 
values \x ± v F t\, keeping in mind that t is slightly shifted to the lower half-plane. As we will 
see very soon, such form is more convenient for taking the Fourier transform. 

The space and time Fourier transform of the above correlation function in the vicinity of 
the one-particle excitation threshold, i.e. at the point (k,u) = (k p ,e p + 5u) (see (|4.4|) ). is 
therefore given by the following quantity: 

dx dt e^-^ { 0\x, t) O(0, 0) ) kp , £p+S u, = p(X) | 2 

e ixup'(X)+it[8ui~ue'{X)] 



x du dxdt — — -, (4.9) 

[-2m{x-v F t)]^-^[2m(x + v F t)]^.-^ 
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in which we have set J\ = [A — es u , A + e$ u ], for some > such that both es^ and Su/esui 
tend to zero in the limit 5lo — > 0. 

Similarly, in the vicinity of the one-hole excitation threshold, i.e. for (k,u)) = (kh,Eh + Su) 
given by (|4.5p . we get 



/ 



dx dt e^*-^) ( 0\x, t) 0(0, 0) ) kh ,e h+ 5u = p(A) 



O |2 



ia;iy(A)+it[5w+i/e'(A)] 



x / dv I dxdt r ; z — -. (4.10) 

[-2m(x - v F t)]( F -^ +1) [2m(x + v F t)]( F -^ 



Integrals over x and t in (|4.9p and (|4.10p are therefore of the type 

„— iPx+iEt 



I a +a (E,P)= / dxdt- = — (4.11) 

+ ' V ' ' J [-2m(x - v F t + i0)] a + [2m(x + v F t - i0)] a - K ' 

- H(E + v F P) H(E - v F P) [47rUjp]a++a __i r(Q+)r(Q _ ) ' ( 4 - 12 ) 

where H denotes the Heaviside step function. The expression ()4.12p for ()4.1ip can easily 
be obtained through the change of variables (<p,ip) = (x — v F t,x + v F t) which separates the 
two-fold integral (|4.11|) into a product of two integrals: 

with 

/ e -ibip r e ibip 

d ^r^FTw = J d1p [2«w-io)]° - (4 - 14) 

R R 

The latter can be easily calculated via substitution 

1 f'°° 

[Tity ± i0)]~ a = — - / e^'a "" 1 da, (4.15) 
"(a) Jo 

leading to 

I a (b) = H(b) . (4.16) 

4.1.1 Behavior of the Fourier transform of the two-point function close to the 
particle excitation threshold 

To obtain the behavior of the Fourier transform in the vicinity of the one-particle excitation 
threshold, we substitute in (|4.12|) : a+ = (i 7 ' 1 4 "_ 1 — l) 2 , a_ = (i 7 ^ -J 2 , and E = 5u — ue'(X), 
P = —up'{\). Hence we have 

E±v F P = 5uj-ve'{X)Tv F vp'{\) = 5u - up'(X)[v ± v F ] . (4.17) 
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Above, we have used the expression (|3,3p for the sound velocity v = v(\). The two Heaviside 
functions restrict the domain of integration on v such that 



v < 



8u) 



p'(X)[v ± v f 



i.e. 



v < 



v < 



5u 



p'(\)[v + v F ] 
Suj 



if 5u > 0, 
if 5oj < 0. 



(4.18) 



p'(\)[v - v F ] 

In (|4.18p . we have used the fact that v > v F (recall that A > q). Hence, 

2vrp(A)|^ 1 ° 1 | 2 



J dx dt e^"**) ( O^x, t) 0(0, 0) ) kp , £p+ 8u 



[4TTv F ] a ++ a — 1 r(a+) r(a_) 



p'<A)[«— tip] 



x [<5w - vj/(\)(v + v F )] a +- 1 [Su - vp'(\)(v - w F )] a -~ 1 di/. (4.19) 

Changing variables u; = (u — t> F )[<5o; — vp'(X)(v + u^,)] /(2v F 8uS) in the first integral, and 
u> = — (f + w F )[5u; — up'(X)(v — v F )]/(2v F Su) in the second integral, we obtain 



O I 2 I Suj \ a++a--l 



Ep + &UJ 



1,-1 I I 2tt ! 



r(a+)r(a_)(v-u F ) a + (v + ^ F ) c 



J dxdte^ t - kx \o\x,t)O(0,0)) kp 

, +oo +oo v 

x\h{5lu) J w a +- 1 {l + w) a -- 1 dw + H(-5uj) J w a --\l + wf+^dwl, (4.20) 

^ ' 

where we have used p'(X) = 2irp(\). Note that we have already taken the Su> — )■ limit inside 
of the remaining integrals, hence extending their integration domain up to oo (recall that es u 
is such that eg u /\6u)\ — > +oo). We stress that, since these limiting integrals are finite, such an 
operation does not affect the leading Suj — > behavior of the overall expression. These last 
two integrals correspond to representations of the beta function: 



w^il + wf^dw 



tt-Hl - -r^du = B(a, l-a-P) = ^ ^ " ° " fi) 



r(i-/3) 



As a consequence, we obtain 

dxdte i{ult - kx \0\x,t) O(0,0)] 



r(i -q+-q_) | jf_J 2 



j.H~(<5u;) sin 7ra_ + <5w) sinvra+j 



7r(u — f F ) a + (v + f F ) a ~ 
a++a- — 1 



(5a; 



2tt 



+ o(|ffw| 



a++a_— 1 



). (4.21) 



- ^2 



We recall that here a + = (F^_ 1 — 1) and a_ = (F 1 _ 1 
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4.1.2 Behavior of the Fourier transform of the two-point function close to the 
hole excitation threshold 

We proceed similarly to obtain the behavior of the Fourier transform (|4.10p in the vicinity of 
the one- hole excitation threshold. We use again (|4.12[> with a+ = (Fj^ l + 1) 2 , a_ = (FTi i) 2 , 
and E = 5u + ve'(\), P = vp'{\). Hence we have 

E±v F P = 5oj + ve'{\)±v F vp'{\) = 5uj + vp (X)[v ± v F \. (4.22) 

Note that, in the present case, v < v F , and the two Heaviside functions restrict the domain 
of integration on v such that 

<u< — — -. (4.23) 



p>(X)(v F +v) p'{\){v F -v) 
In particular, the domain of integration in not empty only if 5lo > 0. Hence 



dx dt e^- k ^ ( 0\x, t) 0(0, 0) ) kh ieh+Su 



2^)1^ | 2 



[Airv F } a + +a -- 1 T(a+) T(a_) 



x H(6u>) [ p ' {x)[Vf v] [Su + vp'{\){v F + v)}^- 1 [Slo - vp'{\){v F - v)} — 1 dv. (4.24) 

/ _ Suj 

Changing variables w = (v F — v)[5u + vp'{\)(v F + v)]/(2v f 5uj) in the integral, we reduce it 
to the integral representation of the beta function, and finally obtain 



dx dt e***-"*) ( Q\x, t) 0(0, 0) } kh ,e h+ 5u, 

= H(Suj)— — ' {27v) . -. r— / w^^Cl-wr-^dw (4.25) 

V T(a+) T(a_) (v F — v) a + (v F + v) a ~ J K ' K ' 

o 

T(a+ + a_) (v F - v) a + (v F + v) a - \2n J 
We recall that here a + = (F^ 11 + l) 2 and a_ = (F~ x 1 ) 2 . 

4.2 Density structure factor 

We now turn to the computation of the density structure factor (|4.ip . As in Section [373^ it is 
convenient to introduce the generating function for the correlation function of densities (see 
formula §32$ ): 



S{k,uj) = -= / dx dt e'M-fc") °° ( e ltH e 2niaQ ^ e~ itH « 

8tt z J ox 1 oa z 

= ^ [ dxdt e *M-M J!_ { e itH e 27riaQ ^ e~ ltH <* ) . (4.27) 
8tt z J oa z a=o 
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Summing over a-twisted form factors of the type (|3.44j) . we obtain 



S(k, u) 



k 2 p(X) 
8vr 2 



d 2 



I fa 12 ix[2ak F +rup' (X)]+it[Su}-Tue' (X)} 
I t, — r I e 



da 2 [-2m(x - v F t)] a +^ [2vri(x + v F t)) a -^ 



, (4.28) 



0=0 



where r = +1 for contributions around the particle threshold and r = — 1 for contributions 
around the hole threshold. The edge exponents read a±(a) = [(a — t)Z + t<1>-|-] 2 , where 
Z is the value of the fractional charge on the Fermi boundary and where $± are given by 
(|3.32p . (|3,33p upon the exchange Ao A. We remind that |J r "_ r | 2 is proportional to a 2 (see 
Appendix [B]). Hence, the second a-derivative has to act on this factor only, meaning that 
we can directly set a = in the remaining part of the integrand. Therefore, re-expressing 
9a[^±i T i] a=0 with the help of (|B.14|) and applying the general scheme described above, we 
get that the leading power-law behavior of the density structure factor on the particle and 
hole thresholds is respectively given by the following contributions: 



S(k, u] 



part 



T(l — a + — a_) Iff _ : 
tt(v — v F ) a + [y + v F ) a 



^H{8uS) sin7ra_ + H{— 5u) sin7ra + | 



Slj 



2tt 



a++a- — 1 



, (4.29) 



around the particle threshold, and 
S{k,u) hoic = H(5u) 



T(a + + a-)(v-v F ) a + (v + v F ) a - \2tt J 



(4.30) 



around the hole threshold. In both expressions (|4.29p and (|4,30p . the edge exponents are 
equal to a± = [Z — <3?±] 2 . They can be easily expressed in terms of the fractional charge and 
dressed phase, 



2 2Z 
Z 1 



2 2Z 

These results agree with the non-linear Luttinger liquid predictions 3^, 3 



(4.31) 
(4.32) 



4.3 Spectral function 

Let us now consider the spectral function (|4.2p . Using the spatial and temporal translation 
invariance of the model, we can recast the initial definition (|4.2|) as: 



A(k,oo) = Slg ^ J dxdt^e i ^ t - kx H^{x,t)^(0,0)} + e' i ^ t - kx \^(x,t)^{0,0))Y (4.33) 



The limiting behavior of each of these two terms can easily be deduced from the general 
scheme described above. 
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The first term (which is non-zero only when oj > 0) gives the leading power-law behavior 
of A(k,u) in the vicinity of (k p ,£ p ) (particle threshold) or of (kh,£h) (hole threshold) as 
explained in Section f4.il We obtain the following result around the one-particle threshold: 



A(k,u] 



part 



(p) (p) 
2^(v-v F ) a + (v + v F ) a - 



oj) sin ira^ + H(—5oj) sin ira^ j 



with exponents 



a 



(p) 

± 



The result is slightly different around the one-hole threshold: 



A(k,u) 



with exponents 



H{Su) I^/J 2 



2vr 



a W +a Cp)_l 



CO i O) , 



(4.34) 



(4.35) 



(4.36) 



(ft) 
a± 



Z± — -$± 
2Z ± 



|±I + #± g> A) 



(4.37) 



The second term of (|4.33p (which is non-zero only for non-positive u) corresponds to 
a Fourier transform of the type (j4.3|) evaluated at the point (— k, — u). Hence, taking into 
account the symmetry k —> —k, the study of this term through the line of Section [4.11 gives 
the leading power-law behavior of A(k,u)) in the vicinity of (k p ,—e p ) or of (feft,— Eh). We 
obtain respectively: 



A(k,u] 



r(l- a f ) -a^ ) )|^_ 1 | 2 

part — (h) (h) 

2ir 2 (v-v F ) a + (v + v F ) a - 



x < H(—oj — e v ) sin ira^ + H(u + e v ) sin ira 



UJ + Er 



for the leading power-law behavior around (k p , —e p ), and 
A (k,^)i^: - 



H(-u;-E h )\^J 2 



2irT(a { l ] + a W ) (v F - v) a + (v F + v)«- 



2tt 



0J + £ h 



(h) , (h) , 



(4.38) 



2?r 



(4.39) 



for the leading power-law behavior around (kh,—Eh)- Note that the edge exponent in ()4.38p 
(resp. in (|4.39p ) coincides with the edge exponent in (|4.36p (resp. in (|4.34j) ) . 

5 Conclusion 



To conclude this paper, it seems appropriate to outline once again the main features of our 
method before discussing possible further developments. In contrast to what happens in the 
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case of massive models, one encounters serious difficulties when studying correlation functions 
from their form factor series in critical models. As we have shown, a macroscopic description 
of the excited states is not enough for a successful summation of the form factor series, and 
a formal replacement of the sum by an integral over the particle/hole rapidities leads to 



meaningless results. This problem was already noted in the literature, see e.g. 78-80]. 

In fact, the success of our method comes from the fact that it is built on a microscopic 
description of the excited states. Such a microscopic description enables us to explicitly sum 
up the quasi-critical form factors belonging to a given class. This amounts to dress the original 
bare form factor (taken for example as the simplest representative of the given class) by the 
complete cloud of excitations having zero momentum and energy in the thermodynamic limit. 
Then, the resulting dressed form factors admit a purely macroscopic description, at least for 
the examples that have been considered above. It shed some light on the relation between 
the particle description in the microscopic model and the effective one corresponding to its 
thermodynamic limit. 

In the present paper, we have dealt with a situation with (contributing) excited states of 
particle/hole type only, and where the particles and holes separated from the Fermi boundaries 
are localized in a small vicinity of a certain point. A possible direction for further developments 
would be to apply this method to more general classes of form factors. In particular, a very 
interesting example to consider is the XXZ Heisenberg chain, for which we also need to take 



into account the contribution of bound states (see e.g. [7l|, |8l|). 
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A Thermodynamic limit of the model 

In this appendix, we remind some definitions of thermodynamic quantities that are useful for 
our study. 

We recall that, in the thermodynamic limit, the set of Bethe roots Xj (|2.2p for the ground 
state densely fills the Fermi zone [—q, q] with a density p(\). This function solves the following 
integral equation: 

P(A) - ^ J K{\ - p) p{p) dp = i-, with K(X) = 9'(X) = (A.l) 
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The dressed momentum p(X) corresponds to the antiderivative of the density function that 
vanishes at the origin: 

q A 
p(A) = po(A) + f 0(X - p) p(p) dp = 2n [ p(p) dp. (A.2) 



-q 

The value of the Fermi boundary q is determined by the dressed energy e(A). The latter is 
defined as the unique solution to the integral equation 

q 

e(A) -^Jk(X-p) e(p) dp = e (X). (A.3) 
-q 

The parameter q is then chosen in such a way that the unique solution to the above equation 
also verifies s(±q) = 0. The solvability of this non-linear problem has been established recently 
in (82|. Note that the value of the parameter q then fixes throught p(A)dA = D the value 
D = lim(A r /L) of the average density. 

Other important thermodynamic quantities are also defined as solutions of linear integral 
equations. Among them, let us mention the dressed phase <j> and dressed charge Z: 

q 

<KA,M)-^ j K{X-u)4>(uJ,p)duj= 9{X ~^\ (A.4) 

-q 
q 

Z(X) - — [ K(X - u) Z(u) du = 1. (A.5) 
2vr J 

-q 

Evidently the dressed phase and the dressed charge are not independent functions. In partic- 
ular, it follows from (|A.4j) . (|A.5|) that 

Z(A) = l + 0(A,-<7)-0(A,<z). (A.6) 

There exists one more useful equation relating Z = Z(±q) and the dressed phase [83| : 

Z- 1 = 1 - <t>{q, -q) - <j){q } q) = 1 + (j){-q, q) + <f>{-q, -q). (A.7) 

We associate to each excited state (defined through the Bethe equations (|2.4p ) a counting 
function It realizes a one-to-one correspondence between the Bethe roots pe- and the 

integers if •) = tj/L. In finite volume, each counting function depends on the particular 
configuration of all integers £j labeling the corresponding state but, in the thermodynamic 
limit (up to corrections in 1/L), the counting functions of any particle/hole states (with N' 
such that N' — N remains finite) coincide. Their common limit is given by the following simple 
combination of the dressed momentum p and average density D: 

lim5H = ^ + ° (A.8) 

L— s-oo Z7T Z 

For a given particle-hole excited state, the counting function allows one to define the particle 
and hole macroscopic rapidities p Pa and ph a as the pre-image of the corresponding integers 
p a and h a - 

ha=CHPa/L), p ha =C\K/L), a = l,...n. (A.9) 
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In general, the existence of particles and holes generates a global shift in the position of 
the Bethe roots of an excited state (|2.4|) with respect to those of the physical ground state 
(E2D: 



n a -\ a = ^^- + 0(L- 2 ), a = l,...,N, a + h x 

Lp{X a ) 



h 



(A.10) 



Here F is the shift function, which depends in fact on the macroscopic rapidities ix Va and fih a 
associated to the positions of the particles and holes: 



F(X) =F A 



0^} 



(A.11) 



It is defined as the solution of the linear integral equation 



AN 



1 



O(X-q) 



TT 



k=l 



2tt 



It can be expressed, in terms of the dressed phase (|A.4j) and dressed charge (|A.5j) . as 

■z{\) 



F(X) = -AN 



(A.13) 



a=l 



with A A?" = N' — N being finite in the thermodynamic limit. 



B Form factor expansion for the correlation function of den- 
sities 



It was shown in |16l that 

e-™ T °*F*m), (B.l) 



(^({M})|e- 2 ™SW|^) _ 



\w*m)\\ w\ 

for an arbitrary state | ipadfi}) ) depending on parameters {/I} satisfying the a-twisted Bethe 
equations f)3.43j) . Here J :a {{'jT\) is given by (|3.45p . and V cx = ~ ^j)- If cu is a real 

parameter, then all ju^, are real, and taking the hermitian conjugation of ()B.2j) we obtain 

\mm)\\m\ " (W) - ( } 

Let us act with the operator d x d a at a = on the l.h.s. of (|B.2jl . Using d x Q[x) = j(x, 0) 
we have 



o x o a ■ 



IK({£})IIIIV> 9 



. (V> g |j(x,0)|^({^})) 



where | ifj'({Jl}) ) = \ ipadfi}) )\ a =o- On * ne other hand, applying the same differential opera- 
tor to the r.h.s. of (|B.2[) we obtain 



nn {^ g \e 2maQ ^w a m)) 

o x o a ■ 



\w«m)\\ u B 



a=0 
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(B.4) 

Here 8^/^ = 1 if {/I} = {A} at a = and 5^,^, = otherwise. We also have used that 
V ex = 2-iraD for the state with £j = j (see (|3.43p ) . Comparing (jB.3[) and (|B.4|) we find that 

(WlMM -^ . ^. ^-V.^ p e ixVc *d F a (\u\)\ (B5) 
IIV>'({m})II ll<M -^'-^^ + 2vr /ex ° UMiJU =0 - l^J 

Similarly the equation (jB.ip leads to the following representation: 

Using these formulas we can prove (|3.42p . On the one hand, we have 
(i(M)i(0,0„= £ e ^ Wj(*.0)l*(<^ 

i*i? }) > ll</-'(M)ll 2 IWI 2 

where the sum is taken over the eigenstates of the Hamiltonian with periodic boundary con- 
ditions, which means that the parameters ju^ satisfy the equations (|3.43p with a = 0. Sub- 
stituting ()B.6p and (|B.5|> into (|B.7p . we obtain 

(i(x,t)i(0,0)) =^ + ^2 E ^ 2 x e-^ + ^(^({/I})| Q=0 ) 2 . (B.8) 

l^({2})>/liM 



On the other hand, one can expand the generating function ( e ltH g^iaQix) e } as 

e-^| <({/!}) 



/ itH 2maQ{x) -itH a \ V- ( ^ | e^^Q(gg-i^ | ^ ^gjj } ( ^ | % ) 

{ ' I II,/,/ (WW 112 ||,/>J|2 ' 



K ({£})) 

where the sum runs over the complete set of twisted eigenstates | )• Using (|B.2p . we 

obtain 

(jtH e 2maQ(x) e -itH a } = ^ ^itE^+ixV^ (^({ju})) 2 . (B.10) 

Note that, in (IBTTOl) . £ ex and P cx are both given in terms of the roots ju^, of the twisted Bethe 
equations (|3.43p . It is easy to see that, acting with the operator —d 2 d^,/8Tr 2 at a = on the 
equation (|B.10p . we obtain ()B.8p . Indeed, 

since the scalar product of two different eigenstates vanishes. Hence, for the corresponding 
terms in the sum (jB.lOp . the second a-derivative acts only on {J ra ({fi})) , while in the rest 
of the formula one can set a = 0. As for the action of the differential operator on the term 
with | ipadfj.}) )| _ = | ip g ), it gives D 2 , since in this case V ex = 2naD. 

Let us now express, for | ip' a ({fl}) )| _ 7^ I )> the derivative ^^"({ju})^^ in terms of 
the form factor of the density operator. Differentiating (|B.2p with respect to a at a = and 
using (|B.lip we obtain 



1) ^ ({/,})| a=Q - 2- pT^^y • (B.12) 
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In particular, if | Y>'({/x}) ) = | if>' rn(h£ ), then 




1) d a P n 



Tn ,t\ a =o 



■a 



= 2m r 



.Q(x) 

rno,£ 



(B.13) 



Taking the derivative over x at x = 0, we finally get 



v a ~r T n Q ,e\ a=0 




rn ,£ 



(B.14) 
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